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Abstract. 

We construct analytical phase-space solutions for pertur- 
bations of flat disks by performing a power series expansion 
for the radius and the velocity coordinates. We show that this 
approach translates into an elegant mathematical formulation 
which is easy to use for a wide variety of distribution functions, 
for as far as resonances do not play a role, such as is the case 
for potentials which are close to quadratic. 
ON As a testcase, the method is applied on the Kalnajs disks. 
ON The results obtained are in full agreement with the analytical 
^~ ^ solutions of the mode analysis. 
Q The strongest advantages of this method are its indepen- 
dence of the mathematical complexity of the unperturbed dis- 
tribution, the degree of detail with which the solutions can be 
calculated and its computational straightforwardness. On the 
contrary, power series solutions are not suitable for describ- 
ing regions where resonant orbits occur, which we therefore 
exclude in this paper. 

We used the technique to analyse perturbations in the cen- 
tral regions of a galaxy, tracking the dynamical consequences 
of a Galactic bar on the kinematics of the solar neighbourhood 
(Hipparcos). We showed how the orientation and strength of 
the bar is related to the properties of the velocity ellipsoid in 
our model. 
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1. Introduction 

The giant spiral structure shown by many galaxies is one of the 
most intriguing features of strongly flattened disks. Through 
the early work of Lindblad, it became clear that these pertur- 
bations were in fact a kind of rotating density waves. Toomre 
(1964) and Lin & Shu (1964) introduced the WKBJ approx- 
imation as a mathematical tool for the description of spi- 
ral structure. They initially supposed the perturbations to be 
quasi stationary and tightly wound. N-body simulations are 
a very valuable alternative, because they provide full, time- 
dependent information on the distribution function. 

On the other hand, mathematical solutions for the per- 
turbed Boltzmann equation can offer additional insights and 
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are sometimes easier to interpret, especially if the mathemat- 
ics are not too dauntingly complex. This equation determines 
the response of a certain unperturbed star distribution (char- 
acterized by a distribution function and a global potential) to a 
perturbative potential. In the self-consistent case, one assumes 
that the perturbative potential and the gravitational potential 
generated by the perturbed mass density are the same, which 
is a legitimate assumption if one is looking at the complete set 
of stars. Unfortunately, it turns out to be almost impossible 
to obtain analytical and general results for this problem (with 
the exception of a few isolated successes, (e.g. Kalnajs, 1971, 
1972), because of the complexity of both equations. In any case, 
Kalnajs (1977) outlined a general method, using action-angle 
variables, with which such analyses could be performed. 

In this paper, we explore a technique for solving the lin- 
earized transport equation for thin disks using a series expan- 
sion form of the perturbed distribution function, which was 
first proposed by H. Dejonghe (1984). The radial coordinate 
and the velocities are expanded in power series, while the an- 
gular coordinate and the time is expanded harmonically. One 
can easily see that this kind of solutions is well suited for the 
central parts (small radius). In this sense, our present approach 
is complementary to the WKBJ approximation, which is de- 
signed for the outer parts of the disk. 

It should be plain from the outset however that power se- 
ries solutions produce regular distribution functions, and have 
a radius of convergence extending only up to the nearest singu- 
larity in the complex plane (the resonances), which are always 
present in any linear approximation. Therefore, we cannot as- 
pire to solve in all its generality the problem of performing a 
mode analysis for a thin stellar disk, but we aim at developing 
a formalism that is easy to implement where applicable. We 
will illustrate the method on such a case which appeared in 
literature (section 3). In addition, we will apply this approach 
on a more realistic distribution function, based on the same 
potential (section 4). 



2. Method 

2.1. Transformation of the transport equation 

Using polar coordinates, we define the distribution function 
for two-dimensional disks as df (r, 6, v r , v' e , t) and the potential 
as V(r,9,t). The corresponding transport equation, written in 
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study perturbations of this equation by sub- 
series expansion, it is convenient to choose 
the expansion as efficient as possible. Circu- 
zero point of choice for the expansion, since 
structure of the resulting equations and most 
galaxies move on more or less circular orbits, 
ular velocity on a distance r from the centre 
equation is obtained by the transformation 
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k = a/4Q 2 + 2rQ^. 

dr 

The symbols a r and ae represent the accelerations caused by 
the perturbing potential, while k of course is the epicyclic fre- 
quency. 

2.2. Power series expansion of the velocities 

Our central assumption concerns the form of the distribution 
function: 
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When we substitute this form in (3) and collect all the terms 
with the same power, we obtain a set of equations that is con- 
veniently represented using vector equations (containing dif- 
ferential operators). To this end we rearrange the coefficients 
in vectors with constant "degree" n = i + j: 
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The transport equation is now equivalent to the set of equa- 
tions (for all n > 0): 
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with the matrix operators defined as (in the following sections, 
matrices are defined by enumerating the nonzero elements) 
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2.3. Harmonic expansion of the perturbation 

In a mode analysis, the time dependence and angular depen- 
dence of the perturbation is considered to be harmonic. The 
potential is thus given by: 



V(r, 9, t) = V m (r)e' (wt - me) + V°(r), 



(11) 



where V° stands for the potential generated by the unper- 
turbed distribution (in the following sections, parameters con- 
cerning the unperturbed situation are written in boldface). A 
more general perturbation can be seen as a superposition of 
such terms, each with different order m and angular velocity 
lo. Because we will linearize the transport equation for small 
perturbations, it is sufficient to examine the behaviour of each 
individual term. Moreover, this linearization and the structure 
of the transport equation allow us to write the corresponding 
perturbed distribution in the same way: 



Pn(r,0,t) = P;(r) e '<" 1 - me) + P°(r). 
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Again, P° is the expansion of the unperturbed distribution, 
while V m and P m are associated with the perturbation. Defi- 
nitions (4) now become 
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and the transport equation (7) transforms into 
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Obviously the unperturbed distribution is already a solution 
of the transport equation, which implies that 
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If we neglect the second order term "D„P n n +1 , a linearized trans- 
port equation for perturbations results: 
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The explicit t and 9 dependence now allows us to calculate the 
effect of the derivation with respect to t and 9 in A„ and B„. 
We define 
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with the new matrix 
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S.^. Power expansion of the radius 

Equation (20) enables us to determine the vectors P™ recur- 
sively, at least in principle. Since 6™ and X>™ contain differen- 
tial operators, it is still difficult to do this numerically. More- 
over, this equation contains a negative power of the radius r, 
which implies that the solution can be divergent in the centre 
of the galaxy, a rather unphysical behaviour. In order to get 
rid of all these problems, we introduce our second assumption, 
which is a power expansion solution for the radius r (of course, 
such a solution cannot contain any singularities in the centre): 
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We call the parameter p the "order". Once more, we define 
new matrices: 
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Since the power expansion of ^ i s given by 
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Collecting the terms in (20) with the same order p, using (30), 
(31) and (33), results in equations of the following form: 
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valid for p = — 1, which we must consider since there are terms 
in Since for m ^ 0, y m (#) = e 8mS is discontinuous in the 
origin, one can assume that a™ = 0. This results in 
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Supposing that the matrices A„ are regular, equation (36) can 
be used for the recursive determination of the vectors P™ p . It is 
convenient to define a new index, the "total order" I = n + p. 
Following increasing total order I and, within a constant I, 
following increasing degree n and decreasing order p (so that 
n + p = I), it is possible to determine all vectors P™ p (see fig. 
1). This holds since, for a particular set n,p with total order 
I, the right hand side of equation (36) only contains vectors 
P™ p i with total order n' + p' lower than I, or with total or- 
der I, but with degree n < n. All these vectors are already 
calculated earlier. One can thus see that (36) enables a unique 
determination of all P™ p . 

However, problems are caused by equation (38), which con- 
stitutes an additional constraint for every total order I, in the 
case n = l,p = 0. Since we saw that (36) already determines 
the solution uniquely, there is no guarantee that the whole set 
of equations (36) and (38) has a solution. In order to obtain 
solutions for this overdetermined system, one can expect that 
constraints on the perturbing potential will appear. In the Ap- 
pendix, we explicitly show that there exist results only when 



fi q = if q is odd, 
a™ = if q <| m |, 



a™ = if q and m do not have the same parity. 



(39) 
(40) 
(41) 




Fig. 1. Scheme for finding the vectors P n ,p recursively. The full 
lines represent the ordering for the determination of the subsequent 
vectors, while the dotted lines indicate which vectors are needed to 
calculate a particular vector. 



These extra constraints are essentially the result of the implicit 
rejection of solutions that diverge in the origin. The first con- 
dition is always fulfilled if V is an even, regular function since 
Q = \J\^r- The second property is similar to the constraints 
on the potentials that are solutions of the Poisson equations 
for a disk distribution, e.g. in (Hunter, 1963). 

Equation (36) now offers a recursive method to obtain the 
vectors P™ p , supposed that the perturbing potential fulfills the 
desired conditions and all matrices A. n (va,s) are regular. This 
regularity requires that vq has a non-integer value (see Ap- 
pendix), which corresponds to non resonance cases. Of course, 
it is not sure that the obtained series is convergent, but one 
can expect that there should be at least some region around 
the circular velocity for which this is the case. 



3. A testcase 

The method presented above essentially solves the collision- 
less Boltzmann equation in the non-resonant regime. In order 
to apply it, one therefore needs an unperturbed distribution 
function, an -ideally- quadratic potential and a perturbing po- 
tential. Such models are not common in literature. Neverthe- 
less, we are in the position to test the method using the Kalnajs 
disks as unperturbed distribution. For these models, the per- 
turbation analysis was already performed fully analytically by 
Kalnajs (Kalnajs, 1972). These models constitute a convenient 
reference point to check the validity of the current calculations, 
since they are uniformly rotating, which avoids problems with 
resonances. In addition, the perturbed potentials are of a poly- 
nomial form and should thus be very well suited to power series 
expansions. 

The distribution function of these disks has the following 
form: 
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fio stands for the angular rotation velocity. Without reproduc- 
ing the mathematics, we can summarize that, by using the in- 
tegral solution for the linearized Boltzmann equation, Kalnajs 
was able to show that for all m, i with i > the potentials 
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create a perturbed distribution of the form (expressed in a 
frame rotating with velocity Q and assuming that fio = 1 and 
d = l) 
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where z = x + iy = re' is a complex notation for the position 
coordinates of the point in which the perturbation is evaluated 
and z = v r + ivg is its derivative. In addition, we have defined: 
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Combining these expressions with Hunter's set of density pairs 
(Hunter, 1963), Kalnajs further proved that the following 
potential-density pairs 
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(with _P ; m the associated Legendre functions and a a constant 
factor) are, for certain values of lo, self-consistent perturbation 
modes for the Kalnajs disks, and thus solutions of the linearized 
transport equation. 

It is not difficult to show that Kalnajs' distribution function 
(45) and our formalism give identical answers. For example, the 
perturbed distribution in the centre is found by setting z = 0. 
This allows to simplify the velocity-dependent part of (45) to 
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where s is a multiplicative factor which is not further defined, 
but in any case independent of z. Since i m = e lm (v r -\-ivg) m , 
it is easy to check that the vectors _P ; ™ (i.e, the coefficients 
of zero radial order p, which describe the distribution in the 
centre) fulfill the condition (38) by comparing (52) to the gen- 
erating function of the eigenvectors of Ai(v, s) (see Appendix). 
In addition, the calculated series expansion in the centre of the 
distribution has exactly the same form as the series expansion 
of the closed form (52). 



is z = ictz. When a is real, this means all stars without ra- 
dial velocity and with angular velocity a. For these stars, the 
distribution function reduces to 
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Again, the calculated series expansion for this fraction of stars 
matches this closed form expression (taking into account that 
Kalnajs' frame is rotating with rotation speed Q, while our 
frame rotates at a speed 1). 

The correspondence between the series expansion and the 
analytical result can be made more clear if the series expan- 
sion is divided by the function F(r,v r ,ve) = ~^ L ( r ' Vr ' Ve )- 



The following formulae give the result for a few low-order per- 
turbations. As predicted by equation (45), the quotient is a 
polynomial with the same total degree as the perturbing po- 
tential (the numerical values of the coefficients are valid for 

LO = 1.5). 
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Integration over the velocity of the series expansion so- 
lution obtained by our method (in the same way as Kalnajs 
did for (45)) yields the series expansion of the perturbed mass 
density. This expansion turns out to be identical to the series 
expansion of the perturbation calculated by Kalnajs' method. 
This shows that, at least for this case, this method is able to 
generate the correct series expansion of the perturbed distri- 
bution. 

4. Application: the velocity ellipsoid in the neighbour- 
hood of a rotating bar. 

If applicable, our method produces an extremely convenient 
representation of the distribution function. Therefore it is easy 
to study the distribution in detail in the presence of a rotating 
bar. The astrophysical interest thereof lies in the possibility of 
observing the stellar velocity ellipsoid in the solar neighbour- 
hood. Several groups have now made a case for the presence of 
a bar at the galactic centre, and in case the perturbative bar ex- 
tends far enough, it may leave its signature in the distribution 
of stars. It is well known that the velocity ellipsoid in equilib- 
rium models does not show any vertex deviation, i.e. the orien- 
tation is always directed towards the centre. The presence of a 
bar may cause a vertex deviation. This is shown schematically 
in fig. 2. The eccentricity can be variable though, depending 
on the amount of energy in radial and tangential motions. 

Since the distribution of the Kalnajs disks is singular at 
the edge, we illustrate the calculation of the velocity ellipsoid 
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with 
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The potential is still a quadratic one, but the velocity dis- 
tribution is now a regular function, decreasing for stars with 
less binding energy. The distribution function and the potential 
are not self-sonsistent, which is reasonable since dfo represents 
only the stars in the central part of the galaxy, while the po- 
tential is generated by the whole disk and even by the dark 
matter which might be present. In addition, it can very well be 
that the sample of stars which is used to determine the velocity 
ellipsoid are only a small subset of the total set of stars in the 
central region of the galaxy, which may moreover be subject to 
observational biases (e.g. a sample of K giants). For that rea- 
son, we do not require the perturbation to be self-consistent, 
but calculate the response of the star sample given by dfo to a 
"generic" barlike potential of the form 
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Fig. 3 represents the response of the star distribution to the 
barlike perturbation, with / = 0.05. For all following results, 
the pattern speed is lo = 2.5 (measured in terms of circular 
rotation speed, which is independent of the radius). In this 
case, the mass density and the potential of the bar have the 
same orientation. This figure shows the velocity dependence 
of the distribution on a rectangular grid in the plane of the 
disk. Apart from the barlike behaviour of the density profile, 
this image shows details of the distribution function, such as 
velocity anisotropy and vertex deviation. 

Fig. 4 shows the dependence of the appropriate parameters 
on the perturbation strength and the position with respect to 
the bar. This position angle is measured from one of the semi- 
major axes of the perturbing bar, in the sense of the rotation 
of the stars. The vertex deviation is defined as the angle be- 
tween the major axis and the direction of the galactic centre, 
without orientation (see fig. 2). The left panel illustrates the 
eccentricity of the velocity ellipse as a function of perturba- 
tion strength, which is determined as the ratio between the 
perturbed and unperturbed potential. Since the unperturbed 
distribution is anisotropic with e = 0.14, this curve starts at 
that value. The right panel draws the vertex deviation as a 
function of the orientation of the bar and the perturbation 
strength. For small perturbations, the eccentricity of the un- 
perturbed distribution is dominating and the vertex deviation 
shows only small deflections from the galactic centre. On the 
other hand, for strong perturbations or almost isotropic distri- 
butions, the vertex deviation becomes an almost linear func- 
tion of the orientation angle. Although the present model is 
not quantitatively correct, it highlights the global dependence 
of the velocity ellipsoid on the bar parameters. 

Looking at the observed stellar velocities in the solar neigh- 
bourhood, there is evidence for the presence of two separate 
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Fig. 3. Representation of the response of the distribution to a barlike perturbing potential. In each point of a rectangular grid on the disk, 
a small density plot of the velocity dependence of the distribution is shown. 




Fig. 4. Left panel: eccentricity of the velocity ellipses depending on the perturbation strength. Right panel: the vertex deviation of the 
velocity ellipses as a function of the orientation of the bar and perturbation strength 
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Fig. 2. Schematic view of the vertex deviation caused by a bar. 
The angle a denotes the orientation angle of the bar, while f3 stands 
for the vertex deviation of the velocity ellipsoid with respect to the 
galactic centre. 



streams (Kapteyn, 1905 and, more recently, Kalnajs, 1991), 
giving rise to a double-peaked velocity distribution. Of course, 
this is incompatible with the simple distribution applied in this 
section. However, since the Boltzmann equation gives is linear 
in the unperturbed distribution, one can simple add two un- 
perturbed components with different streaming velocities. The 
response to the bar potential will show the same components, 
each with its own vertex deviation pointing in the same direc- 
tion. 



5. Discussions 

The series expansion paradigm for solving the linearized trans- 
port equation without resonances yields an elegant mathemat- 
ical description, with an easy computational algorithm as a re- 
sult. The input required to start the calculations is the power 
series expansion of the time-independent and axisymmetric dis- 
tribution, together with the power expansion of the harmon- 
ical perturbing potential. The calculations themselves consist 
of very simple algebra, which can easily be done by a com- 
puter and turns out to be very fast. One thus can see that this 
method is in principle applicable to all distributions, regardless 
the mathematical complexity. This is presumably its strongest 
advantage, since analytical perturbation analyses almost al- 
ways turn out to be very complicated. If this way of solving the 
linearized transport equation is coupled to a method for solving 
the Poisson equation for disks, one should be able to perform 
fully self-consistent mode analyses for more complex distribu- 
tions. Moreover, this method highlights also some mathemat- 
ical symmetries of the distribution, since we can show that 
many terms of the expansion should be zero (see Appendix). 

The fact that the perturbations are represented as a full 
phase-space distribution makes them well suited for detailed 
analysis. As an example, we calculated the vertex deviation 
in the solar neighbourhood due to the presence of a bar. A 
significant vertex deviation should be observable using data 
from the Hipparcos satellite, which therefore will enable us in 
principle to obtain further evidence for the presence of a bar 
in our Galaxy. In addition, the obtained value for the vertex 
deviation could enable one to estimate the orientation of the 
bar. 



become important, since the power series expansion is by no 
means able to handle poles in the solution. Even for a growing 
perturbation (with complex u), the problem still exists since 
the radius of convergence around circular orbits will exclude 
all poles, even in the complex plane. This essentially restricts 
this method to the following cases: 

— Quadratic or almost quadratic potentials. In this case, res- 
onances only occur for particular frequencies and they then 
occur over the whole disk. The solution thus will have 
no poles. For galaxies which are not rotating uniformly, 
a quadratic approximation still makes sense for the central 
region of the system. 

— Outer regions of not uniformly rotating disks, excluding all 
resonant radii. 

— Rapidly rotating or growing instabilities. In these cases, the 
stars suffer no important resonances. 

If our method is applicable, there is a substantial opera- 
tional gain, since computers are extremely efficient in evaluat- 
ing power series. The fact that only an expansion is obtained 
instead of a closed form solution is, in our opinion, a minor 
disadvantage, since often closed form solutions are so complex 
that the only way to handle them is to calculate the series 
expansion! 

It is possible to extend this method to a certain degree of 
non-linearity by using more than one harmonic in the expan- 
sion and calculating all the nonlinear terms in (15) which have 
a harmonical degree lower than m. Moreover, the other nonlin- 
ear terms can be taken into account by a perturbation analysis. 
This we reserve for a future paper. Another important possible 
extension is the introduction of a certain thickness of the disk, 
entering the z axis and the v z velocity in the series expansion. 

A. Appendix 

A.l. Eigenvalues and eigenvectors of A 
We represent the eigenvectors of 



A„(0,s) 



t s 
} 2s 

n — l 



\ 



(Al) 



V i / 

corresponding to the eigenvalue fi by 



( <o( s )' 

X n,l( S ) 



(A2) 



\ x n,n (^) ) 

The elements of the vector thus have to satisfy 
n — i + 1 



-<,_!(*) = 0, i = 0. 



By defining 



xJ 1 ■ (s) = s 1 xJ 1 

n,i V / n,\ 



(A3) 



(A4) 



~Wn,i + (« + 1 ) X n, t + l + {n~i + l)< i8 _! 

i = 0...n. 

The equation for the generating function 

n 

8 = 

follows from (A5): 

rldl 1 rlG^ 

— fi&Z - z — 1- nzG'„ = 0, 

az az 

and has a general solution of the form 



(A5) 



(A6) 



(A7) 



(s + 1)^(1 



Since G has to be of polynomial form, only integer values for 
the eigenvalues fi, with \fi\ < n and with fi the same parity as 
n are valid. It now follows immediately that the more general 
coefficients x^ t are generated by 



G»(s,z) = (- s + 1)^(1 



(A9) 



More generally, one can see that the matrices A„(r, s) have the 
same eigenvectors, but with r + fi as eigenvalues. This proves 
that these matrices are regular for all non-integer values of r. 
Another set of vectors, frequently used in the text and closely 
related to these eigenvectors are 



\ 



(A10) 



which is essentially E k (s) padded with n + I — k zeroes. The 
following identities, which can easily been proven, 

Gt±l(s,z)-Gt(s,z) = ± Z -Gt(s,z) (All) 

(k-zj ; )Gt(s,z) = Tll Gtt\(s,z) 

+ (k±v)G£_ 2 (s,z) (A12) 
translate to corresponding properties between eigenvectors: 



I 







s \K,k( s ) 
(k 

k - 1 



E' 1 

^n+l,k 





(A13) 



+(k±n)E" n 



2 « 



(A14) 



A. 2. Existence of solutions 
A. 2.1. Rewriting the equations 

In this section, we will prove that (36) and (38) do have solu- 
tions under the following conditions for the potential: 

(1) Q q = if q is odd 

(2) a™ = if q <| to | (A15) 

(3) a™ = if q and to do not have the same parity. 



to rewrite the equations in another format. We rewrite all vec- 
tors in terms of the eigenvectors of A n (0, i) (we will represent 



P" 



' ^ ^ ^ ^ C n,p,k ^n,k ' 
k = fi=—k;2 
n k 

n 'P / y / y ^n,p,k ^n,k ' 

k = u=—k;2 



(A16) 
(A17) 



_ — \k 

where the sum sign 2^ M __j, 2 indicates that the summation 
index fi is, starting with —k, each time increased by 2. In the 
following equations, we will denote the double summation over 
(A8) k and /t by ~^2 k 2 Combining the properties derived in the 



previous section, the transport equation transforms into: 
E <p," [[" - "o + (« - - (n - 

k,u;2 



"EE "oC-a,* + f ) " (» " fc )( 2 + 

fc,/x;2 a = 2 

+(« - *)( 2 + «)£„t +1 + «t-^«+ i+£2 - ^ty 

+ (fc + e 2 ^)«+ 1 - 1 -<,-2)]] 

+ E C n-l,p+l,k [i m + ») E n,k 
k,;j.;2 



+Cl (n-p-k-2)(E££ 1 -EZ ik )] 
p+i 

E E C °+l,p-c + l,* a « [( m - - t l ) E n+l, 



+m(n + i _ k)E» +l k + (mti + a*)££Ti,*-iJ , (A18) 

where t\ and C2 are sign constants which should be +1 or — I. 
The constraint (38) now simply becomes 



c To, % = if (/t / -to) or (A; / I). 



A. 2. 2. The unperturbed distribution 



(A19) 



First, we need some general properties of the series expansion 
for the underlying unperturbed distribution, implied by equa- 
tion (16). After rewriting this equation in the same way as the 
perturbed case, this reads: 

v 

2n X(0,i)P°, p + ^S),[2X(0,i)- g iqP°„, p -, 

q=2 

+ B n , P+ iP°_i, P +i =0, (A20) 
and the boundary condition becomes 

^i(0,i)P° o =0. (A2I) 

Though this equation looks similar to the general case, it does 
not allow to determine all P„ iP , because the matrix A is singu- 
lar and there remain some degrees of freedom. But it is useful 
to extract some symmetry properties from it. Using the pre- 
viously derived identities, it is again possible to rewrite this 



ra - k ( 



2 + a / 



/t / \ k + fi / 



J k-2 I 



ra - p - k - 2 / 



(A22) 



In the same way, the boundary condition translates to 



if (/t / 0) or (k / I). 



We will now prove by induction that 



,0,11 



if /t / 0, 



(A23) 



(A24) 



in other words, the unperturbed solution can be written solely 
in terms of the eigenvectors of A n (0,i) with eigenvalue 0. 
Equation (A23) shows that this holds for all ra, with p = 0. 
Supposing it is true for all n',p' with ra' + p < ra + p or 
ra'+p' = ra + p and p' < p, we have to prove that (A24) 
holds also for ra — l,p + 1. Looking at the coefficients of E^ k , 
which is independent of all other vectors in the equation (since 
it is the only vector with a non-zero first element) and using 
the induction hypothesis, we have that c°'^j p+1 k .(i = 0. This 
completes the proof. The equation thus becomes (the summa- 
tion sign indicates that k is incremented by 2): 



o,o n - k ( 

El 



q -<°p-c, J(™ - k ) l v E k 

k;2 a = 2 

Eo,o n — p — k — 2 f 

C n-l,p+l,k 7 



2 + a / \ k I 
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(A25) 



Another property, proved by induction, 



~n,p,k 



if ra + p is odd. 



(A26) 



For all ra, this is obvious for p = because of the condition 
(A23) and since k should always be even (p = 0). Keeping in 
mind the conditions (A15), one can see, using (A25), that this 
holds also for p / 0. Finally, we want to prove that 



~n,p,k 



if k < n — p. 



(A27) 



For all n, with p = 0, (A23) again ensures that this is true. If it 
is true for all pairs (n,p) with ra'+p' < ra+p or ra'+p' = ra+p 
and n < n, the left hand side of (A25) only contains vectors 
for which k > n — p. One can see that this allows to conclude 
that (A27) holds also for n,p. 



We will now first show by induction that 



^n,p,k 



if n + p + m is odd. 



(A28) 



For n = 0, p = 0, this is obviously fulfilled. For all n, p, we 
have to prove that it is true, given that it holds for all pairs 
ra',p', with n +p' < n +p or n +p' = n +p and n < n. Since 
we have (A15) and using the properties of the unperturbed 
situation, one can verify by inspection that the right hand side 
of equation (A18) only deliver vectors E 1 ^ k if ra+p + m is even, 
supposed that £2 is choosen I for fi ^ k and —I for fi = k. This 
is sufficient to conclude that the required property holds for 
the terms n,p on the left hand side. 
The last property which we need is 



^n,p,k 



if fi < —m or k — fi < n — p + m. 



(A29) 



Again, it is easy to check that this holds for n = 0, p = 0. 
As induction hypothesis, we assume that it holds for all n' , p , 
having ra'+p' < ra +p or ra'+p' = ra + p and ra' < ra. Choosing 
t\ = +1 for terms which have k — fi > ra— p + m — 2, t\ = — I in 
the case k — fi = n— p + m — 2 (other situations are impossible 
because of the induction hypothesis and the previous parity 
properties) and using the properties about the unperturbed 
situation, one can see by inspection that the right hand side of 
equation (A18) only delivers terms E^° kg with fio > — m and 
ko — /to >ra— p + m (the same choice as for the previous case 
holds for €2). Each of these terms is compensated on the left 
hand side by the following contribution: 



(k - fi / k - Ho) V (/t < /t ) 
I 



n,p,k 







r m,fi 
~n,p ; kQ 



)io — vo — k — ra 



,,"1,(10 + 1 



k — ra 



(po — vo + k — n)(fio — vo + &o — ra + 2) 



(A30) 
(A3f) 

(A32) 
(A33) 



r m,fi + n — k 
^n,p,n 



(k — n)(k — ra + I ) 



(/t — vo + k — n)(fi — vo + k — ra + 2) 
(-1) 



(/t — vo + ra — ko) 



(A34) 



This shows that the induction hypothesis is also satisfied for 
ra, p. Applying (A29) on p = 0, ra = I, one can conclude that 
c^oX onl y if A* > — m and k — fi > I + m. Combining 
this inequality with the properties /t < A; and k < I, we can 
conclude that c™ '^ ^ only if /t = —ra and k = I. 
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